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ABSTRACT 

Double-peaked oxygen lines in the nebular spectra of two peculiar Type Ib/c Supernovae (SN Ib/c) 
have been interpreted as off-axis views of a GRB-jet or unipolar blob ejections. Here we present 
late-time spectra of 8 SN lib, lb and Ic and show that this phenomenon is common and should not be 
so firmly linked to extraordinary explosion physics. The line profiles are most likely caused by ejecta 
expanding with a torus- or disk-like geometry. Double-peaked oxygen profiles are not necessarily the 
indicator of a mis-directed GRB jet. 

Subject headings: gamma-ray burst: general supernovae: general 



1. INTRODUCTION 

To understand the connection between long-duration 
Gamma Ray Bursts (GRBs) and the supernovae (SNe) 
associated with them, it would be significant to detect 
the effects of jets that are produced in the explosion 
even when the observer is not in the relativistic beam 
(jWooslev fc BloomI 120061 ) )■ In the broad-hned Type Ic, 
SN 2003jd, an off-axis GRB jet was proposed as the 
cause for the double- peaked oxygeii profi l e observed in 
its ne bular spectrum (iMazzali et al.l l2005l : IValenti et al.l 
l2008l ). A similar doulDle-peaked oxygen profile, with a 
large blues hift, was detected i n the peculiar Type lb, 
SN 2005bf (|Maeda et al.l [20071 : iModi ad [200l . where it 
was interprete d as coming from a unipolar blob or jet 
(|Maeda et al.l [20071. 

If the double-peaked oxygen shape were exclusively 
caused by highly relativistic jets, then we might ex- 
pect only broad-lined SN Ic, whose line widths ap- 
proach 30,000 km s~^, and which are the only subtype 
of stripped-envelope SN seen at the sites of GRBs to ex- 
hibit them. However, if we see double-peaked profiles 
in a range of supernova types, we might conclude that 
asphericities are present in more typical core-collapse 
events. 

We set out to obtain a set of nebular spectra for su- 
pernovae of various subtypes to measure the emission 
line profiles with adequate resolution and signal to detect 
signs of these double-peaked profiles in nebular spectra 
of supernovae. As the supernova turns optically thin 
a few months after maximum light, the emission line 
shapes can provide i nformation on the ve l ocity distri- 
bution of the ejecta (iMatheson et"aLl 120011 : iFolev et all 
120031 : iMaeda et al.ll2006D . Structu re in emission line pro- 
files have been reported in SN II (ISp yromili^'igQl*). SN 
IlbJSpvromilio 1994; Mat heson et al . 2000) and SN lb 
(|Sollerman et al.lll998l : lElmhamdi et al.ll2004ft . and inter- 
preted as clumping of oxygen in the ejecta. 
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In § [2] we describe the observational sample. In § [3] 
we discuss the line profiles and show that the double- 
peaked line profiles are unlikely to be caused by optical 
depth effects. We speculate on the implications of our 
results in § m and summarize in § [5l 

2. OBSERVATIONS AND ANALYSIS 

Spectra were obtained with the 6.5 m Clay Telescope of 
the Magellan Observatory located at Las Campanas Ob- 
servatory, the 8.1 m Gemini-North telescope via queue- 
scheduled observations (GN-2005B-Q-11, GN-2006B-Q- 
16 PI: Modjaz), the 6.5 m Multiple Mirror Telescope 
(MMT) and the 1.5 m Tilhnghast telescope at the Fred 
Lawrence Whipple Observatory (FLWO). The spectro- 
graphs utilized were the LDSS-3 (Mulcha ey fc Glad- 
ders 2005) at Magellan, the CMO S-North jH ook et alj 
120031 ) at Gemini, the Blue Channel (ISchmidt et al. 198^ 
at the MMT, and FAST (Fabri cant et al.lll998^ at the 
FLWO 1.5 m telescope. All optical spectra were reduced 
and calibrated employing standard techniques in IRAF^ 
and our own IDL rou tines for flux calibration (see e.g., 
IMatheson et al.ll2008l ). 

In Figure [1] we present the observed nebular spec- 
tra of 8 new SN lib, SN lb and SN Ic. For clarity we 
show only the latest spectrum of each SN (except for 
SN 2004ao). We double the number of unambiguous SN 
lb with late-time data from four historical cases (SNe 
1983N, 1984L, 19901, and 1996N) to a total of eight 
(now including SNe 2004ao, 20 04dk, 2004gq, 200 4gv), 
not counting SN 2005bf (see also lMaeda et al.ll2008f) . We 
determined the p hase with resp ect to maximum from our 
own photometry (|Modiajl2007l . M. Modjaz, in prep.). In 
addition, we show our data of the peculiar SN lb 2 005bf, 
indepe ndently observed and analyzed by .Maeda et al.l 
(|2f)0% . and of SN 2006aj the broad-lined S N Ic con- 
nected with GRB 060218 (iModiaz et al.ll2006D. indeoen- 
dently observed and analvzed bv IMazzali et all (|2007D . 

All spectra display the hallmark features of nebu- 
lar SN Ib/c spectra: strong forbidden emission lines of 
intermediate-mass elements such as [O I] AA6300,6363 
and [Ca II] AA7291,7324, similar to spectra of SN II. 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 1. — Selection of nebular spectra of SN lib, SN lb 
and SN Ic in their respective rest frames. SN name, type and 
phase of spectrum (with respect to maximum light, except for 
GRB 060218/SN 2006aj which is referenced to time of GRB burst) 
are marked. For clarity we show only the latest spectrum of each 
SN (except for SN 2004ao), even though we obtained several for 
some SN. Also the main nebular emission lines of [Mg I A4571, 
[O I] AA6300,6363, [Ca II] AA7291,7324 and O I A7774 are marked 
at the very top. 

In the absence of hydrogen, oxygen is the primary 
coolant in the ejecta of stripped-envelope SN at late 
epochs when the gas is neutral or at most singly ionized 
(jUomoto fc Kirslmei| [T986: Fransson & Chevalier 1987). 
SN 2005bf is the only SN discu ssed here that exhibits 
broad Ha (see also lMaeda et al.li2007) . 

3. DETECTION OF DOUBLE-PEAKED OXYGEN 

LINES 

In the optically-thin case, the late-time emission line 
profile is dictated, in principle, by the geomet ry and dis- 
tribu t ion of the emitting materia l (jWansson fc Chevalieil 
Il987t ISchlegel fc Kirshne^ 11989'). A radially expanding 
spherical shell of gas produces a square-topped profile, 
while a filled uniform sphere produces a parabolic pro- 
file. In contrast, a cylindrical ring, or torus, that expands 
in the equatorial plane gives rise to a "double-horned" 
profile as there is very little low-velocity emission in the 
system, while the bulk of the emitting gas is located at 
±Ut, where vt is the projected expansion velocity at the 
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Fig. 2. — Montage of SN with double-peaked oxygen profiles in 
velocity space. SN name, type and phase of spectrum (with respect 
to rnaximum light). SN 1984L , a SN lb is from Schlcgcl & Kirshnei] 
( 198 9). SN 19941 (Filip penko et al. 1995) that exhibits a simple 
parabolic oxygen line profile is plotted for comparison at the bot- 
tom. The dashed line marks zero velocity with respect to 6300 A 
. For SN 2004ao, we plot the scaled profiles of O I A7774 (in blue) 
and [Mg I A4571 (in red), which are not doublets, but also exhibit 
the two peaks. As discussed in the text (fj 13.11 1. the two horns are 
unlikely to be due to the doublet nature of [O I] AA6300,6363. Note 
the large blueshift in SN 2005bf compared to the other SN which 
we do not explain here. 



torus. 

Figure [2 shows that SNe 2004ao, 2004gt, 2006T show 
conspicuous double-peaked lines of [O I] AA6300,6363. 
We plot the oxygen line profiles in velocity space for 
these objects in Figure O A literature search of pub- 
lished SN lib, lb, and Ic spectra reveals that SN lb 
1984L, the prototype of the SN lb class, also shows a 
clear double-peaked profile that went unr emarked in the 
original publication (jSchlegel fc Kirshnei| [l989l . We plot 
the spectrum with the highest S /N in Figure [S] Finally, 
we also include our data on SN 2005bf, the peculiar SN 
lb with extraordinary early-time light curv es and spectra 
(|Tominaga et al.ll2005l : [Foiatelli et al.ll2006f ). that showed 
highly blue -shifted (by ~ 2000 km s~M oxy gen and cal- 
cium lines CMaeda et al1l2007t lModiadl2007[ ). 

This double-peaked line profile is visible in the oxygen 
line and not in calcium. For SN 2004ao, the only SN 
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for which we have a sufficiently high S/N spectra, the 
double-horned feature is present in the permitted line 
of oxygen, O I A7774, and in [Mg I A4571 (see Fig. H) 
Our late-time spectral follow-up for some SN (iModiad 
HQOT") shows that SN without double-peaked profiles (SNe 
2004gk, 2004gq) do not exhibit them at any epoch, while 
SN with double peaks (SNe 2004ao, 2004gt) retain them 
over the full observed period. 

The two horns are roughly symmetrically offset by 
-f 000-2000 km s"^ around the trough. For most SN 
in our sample, the trough between the two "horns" is lo- 
cated at nearly zero velocity for [O I] 6300A. SN 2004ao 
shows an apparent blueshift and SN 1984L shows a red- 
shift in the trough. A velocity shift of ± 200 km s~^ (as 
seen for SN 2004gt) could be due to rotation of the host 
galaxy. We discuss other possible causes for these shifts 
in §12 

3.1. Optical Depth Effects causing the Double Peaks? 

Could the observed double-peaked profile be due to 
optical depth effects? The [O I] AA6300,6363 doublet has 
a velocity separation of 3000 km s~^ (with respect to 6300 
A) and an intensity ratio of 3:1 in the optically thin limit. 
Optical depth effects for [O I] AA6300,6363 hav e been 
discus s ed in the su p ernov a context by Lcibundgut et al] 
11991), Sp vromiliol (1991) and iLiXMcCravl (11992'). In 
these models, the ratio of [O I] A6300 to [O I] A6363 
evolves from ~1 to ~3 from early times to late times 
as the supernova expands and the lines become op tically 
thin as was observed for SNe 1987A fSpyromili o et al.l 
[19911 : iLi fc McCravl[l99l and 1988A l|S^vromilidlT99li r 
In Figure [31 we plot the measured ratio of the blue to 
the red horn for SN 2004ao, the supernova in our sample 
with the longest time span of observations. The ratio 
of line intensity decreases over time for SN 2004ao. In 
contrast, this ratio would increase with time if the two 
horns were d ue to [O I] A6300 an d [O I] A6363, as seen 
in SN 1987A (|Li fc McCravl[l99l overplotted in Fig. [3]). 

In SN 2004ao, where we have sufficiently high S/N- 
data, we also see a similar double-peaked profile in the 
permitted O I A7774 line and in [Mg I A4571 which are 
not doublets, as shown in Fig. [2] This suggests that the 
line shapes really do provide a guide to the distribution 
of oxygen and magnesium in the ejecta. Optical depth is 
not the probable cause for seeing two equal peaks in the 
[O I] AA6300,6363 line profile. 

4. INTERPRETATION 

We suggest a torus-like structure for either the oxygen 
distribution or for the radio active energy source that ex - 
cites the oxygen lines, s^Co ([Fransson fc Chevaliedll989f ). 
If the distribution of the energy source were most im- 
portant, we would expect to see a double-peaked profile 
for all the emission lines, including calcium. Since the 
double-peaked lines are seen only for oxygen and magne- 
sium, the evidence favors the view that we are seeing the 
effects of the distribution of those elements in velocity 
space. 

Models of the line shape that r esults from an expand - 
ing disk are worked out b y, e.g., iGerardv et al.l ()2000( ). 
and Fransson et al.l ()2005f ). Comprehensive modeling of 
the line profiles needs to take into account viewing angle 
effects and the thickness of the line-emitting and expand- 
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Fig. 3. — Evolution of the intensity ratio of the blue to the 
red horn ( I(Blue Horn)/I(Red Horn) ) as a function of time for 
SN 2004ao (filled circles, this work). The intensity ratio decreases 
over time. In contrast, if the blue and red horns corresponded to 
[O I] A6300 and [O I] A6363 the intensity ratio woul d increase over 
time, as observed in SN 1987A (em pty circles, from IlT &: McCravl 
[1993) and 1988A IISpvromiliol(T99llV 



ing region. We leave that as a topic for future work. 

Here we do not attempt to explain the extreme 
blueshift (- 2000 km s"^) seen in SN 2005bf at epochs 
later than 200 days. The less extreme blueshifts seen 
in spectra obtained at t < 200 days for one of our SN 
(SN 2004ao) might be caused by the same mechanism 
that is responsible for the observed blueshifts in SN that 
show a single-peak line of [O I] AA6300,6363. S. Tauben- 
berger at al. (in preparation) observe in a large set of 
late-time spectra of 34 SN lb and SN Ic that the oxy- 
gen line centroids are found to be blueshifted for spec- 
tra taken at t < 200 days. These blueshifts range up 
to '^1500 km s'^ for spectra at t ~90 days and go to 
zero with increasing time. Taubenberger et al. exclude 
dust formation, contamination from oth er lines and 
geometric effects as potential causes and invoke resid- 
ual opacity effects as the most likely reason. The ob- 
served blueshift-values in our sample agree with those 
seen in Taubenberger et al. for similar epochs; thus, the 
same mechanism might be causing the blueshifts in both 
single- and double-peaked [O I] AA6300,6363. We en- 
courage future detailed modeling of the radiative trans- 
fer of [O I] AA6300,6363 to elucidate the exact reason 
for the blueshift. We do not have a simple explanation 
for the redshift seen in the trough of the double-peaked 
[O I] AA6300,6363 of SN 1984L. 

We conclude that the most probable explanation for 
the double-horned oxygen profiles is a torus like, or at 
least flattened, distribution of oxygen, that may be a 
relic of the explosion physics. Energy deposition and ra- 
diation tra nsport modeling of late-time emission in SN lb 
(jFransson fc Chevalier 1989) predicts that material con- 
tributing to [O I] AA6300,6363 (and to [Mg I A4571) 
emission is situated further out in the ejecta than ma- 
terial emitting in [Ca II] AA7291,7324 (see their Fig. 
8). We speculate that the mechanism causing the as- 
phericities has to affect the outer layers of the SN debris 
more than the inner ones. The anisotropics we observe 
in the oxygen lines suggest large-scale plumes of mix- 
ing rather than the small-scale inhomogeneities invoked 
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to explain sub-structure of the oxygen l ines a s reported 
for SNe II 1979C, 1980K (iFesen et all [l999l and refer- 
ences therein), in SN l ib 1993 J ("Matheson et al.ll2000D 
and in SN Ib/c 1985F (iFilippenk o & SargcnlJ|l98l~In 
the case of our observed SN, the anisotropics may be of 
global nature in order to give rise to such a clear double- 
peaked line profiles. Alternatively, in case the super- 
nova progenitors are part of binaries, binary interaction 
or merger migh t be modulating the geometry of the su- 
pernova ejecta (jMorris fc Podsiadlowskill2007[ ). 

Such a clear signature of a pure double-peaked 
oxyge n profile has not be en found in other SN 
Ib/c (iMatheson et all 1200 If ), except in SN 2003jd 
dMazzali et al. 200Cand in SN 2005bf. SN 20G3jd was a 
broad-lined SN Ic without an observed GRB and showed 
the s ame double-peaked oxygen profile (Mazz ali et al.l 
|2005[ ) as SN 2004ao (c o mpar e their Fig. 2 with our 
Fig. [2]). iMazzali et al.l (|2005f) interpret their observa- 
tions as indicating an aspherical axisymmetric explosion 
viewed from near the equatorial plane, and suggest that 
this asphericity was caused by an off-axis GRB jet. They 
conclude that only special SN, those with high expan- 
sion velocities and possibly connected with GRBs, are 
aspherical, and use SN 2003jd as a link between the nor- 
mal, spherical SN Ic and those highly aspherical ones 
connected with GRBs. In our sample, however, the SN 
showing these line profiles are normal SN lb and even SN 
lib, i.e. SN from stars with intact helium and partial H 
envelopes before explosion, and have normal early-time 
expansion velocities. Thus, it appears that asphericities 
are prevalent in normal core-collapse events. 

Indeed, complex and ring-like velocity structures of 
oxygen have been observed in SN remnants (SNRs), most 
prominently in the oxygen-rich SNR 1E0102. 2-7219 i n 
the Small Mage llanic Cloud (|Tuohv fc Dopital I1983D . 
iTuohv fc Dopital (jT983) and subsequent papers fit a 
twisted ring model to the spatial and velocity extent of 
the filaments that extents to velocities ranging between 
—2500 to + 4000 km s^^. It is conceivable that the stel- 
lar death leading to SNR fE0102. 2-7219 had a geometry 
for the ejecta that was similar to the events we observe. 
Asphericities may be com mon in in core-colla pse events, 
be they neutrino-dr iven (jScheck et al.l I2006D. acoustic 
(iBurrows et al.ll2006D. or magneto driven ([Burrows et al.l 
l2007t iDessart et al.ll2008[ ). as indicated by polarization, 
neutron star kick velociti es and the morpho logies of 
young remnants. Recently, iMaeda et al.l (|2008f l also re- 
ported that double-peaked line profiles are not rare in 
stripped core-collapse SN. They suggested that the ob- 
served fraction of double-peaked profiles (~ 40 ± 10%) 
is consistent with the hypothesis that all of these core- 
collapse events are mildly aspherical. Although their 
infe rences are based on the results of bipo lar jet mod- 
els (|Maeda et al.ll2006l ) . iMaeda et al.l (|2008l ) suggest that 
a wide variety of geometries, similar to those suggested 
here, may be present. 

5. CONCLUSIONS 
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In summary, we detect clear double-peaked line pro- 
files of [O I] AA6300,6363 in three SN lib and SN Ib/c 
out of our observed sample of 8 and additionally, in one 
out of four published SN lb. Our sample does not in- 
clude broad-lined SN Ic or peculiar objects. These line 
profiles are probably not caused by optical depth effects, 
and suggest global anisotropics in the ejecta. Prior to 
this work, double-peaked oxygen lines had only been re- 
ported in two peculiar SN Ib/c and interpreted as off-axis 
GRB-jet or unipolar blob ejections. It seems more likely 
that asphericities are present in a wide variety of core- 
collapse events and they are not strictly confined to the 
supernovae associated with GRB jets. Although special 
models have been proposed to account for the line pro- 
files in peculiar supernovae, our investigation suggests 
that double-peaked profiles, and underlying disks, are 
not unusual. These ordinary SN from our sample have a 
variety of subtypes, reflecting the diversity of mass loss 
prior to the explosion. 

We recommend a meta-analysis of available nebular 
spectra of all core-collapse SN, spanning the full range 
from SN II to lib, lb, Ic and finally, to broad-lined SN 
Ic, in order to quantify the kinematics and geometry of 
the ejecta and to seek trends as a function of SN type. 
Further high-SN/N multi-epoch observations of a sample 
of SN Il/IIb/Ib/Ic coupled with radiative transfer models 
should help to elucidate the observed blue- and redshifts 
of the line profiles. 



M. M. would like to especially thank R. Fesen for 
very helpful suggestions and for reading a draft of this 
manuscript, and in addition C. Fransson, R. Narayan, D. 
Sasselov, B. Leibundgut and J. Spyromilio for insight- 
ful discussions. We thank and FLWO 1.5m observers 
for obtaining service-spectroscopy, and the referee, Alex 
Conley, for constructive comments. M. M. acknowledges 
support from the Miller Institute for Basic Research dur- 
ing the time in which part of this work was completed. 
Observations reported here were obtained at the MMT 
Observatory, a joint facility of the Smithsonian Institu- 
tion and the University of Arizona, at the 6.5 meter Mag- 
ellan Telescopes located at Las Campanas Observatory, 
Chile, at the F.L Whipple Observatory, which is oper- 
ated by the Smithsonian Astrophysical Observatory, and 
at the Gemini Observatory, which is operated by the As- 
sociation of Universities for Research in Astronomy, Inc., 
under a cooperative agreement with the NSF on behalf of 
the Gemini partnership. Supernova research at Harvard 
University has been supported in part by the National 
Science Foundation grant AST06-06772 and R.P.K. in 
part by the Kavli Institute for Theoretical Physics NSF 
grant PHY99-07949. 

Facilities: MMT (Blue Channel spectrograph), 
FLWO: 1.5m (FAST), Magellan:Baade(LDSS3), Gem- 
ini:Gillett (GMOS-N) 



Dessart, L., Burrows, A., Livne, E., & Ott, C. D. 2008, ApJ, 673, 
L43 



5 



Elmhamdi, A., Danzigcr, I. J., CappcUaro, E., Delia Vallc, M., 
Gouiffes, C, Phillips, M. M., & Turatto, M. 2004, A&A, 426, 
963 

Fabricant, D., Chcimcts, P., Caldwell, N., & Geary, J. 1998, PASP, 
110, 79 

Fescn, R. A., et al. 1999, AJ, 117, 725 

Filippenko, A. V. & Sargent, W. L. W. 1986, AJ, 91, 691 

Filippenko, A. V. et al. 1995, ApJ, 450, Lll 

Folatelli, G., et al. 2006, ApJ, 641, 1039 

Foley, R. J., et al. 2003, PASP, 115, 1220 

Fransson, C., et al. 2005, ApJ, 622, 991 

Fransson, C. & Chevalier, R. A. 1987, ApJ, 322, L15 

— . 1989, ApJ, 343, 323 

Gerardy, C. L., Fesen, R. A., Hoflich, P., & Wheeler, J. C. 2000, 
AJ, 119, 2968 

Hook, I., ct al. 2003, in Proceedings of the SPIE., cd. M. lyc & 

A. F. M. Moorwood, Vol. 4841, 1645-1656 
Leibundgut, B., Kirshner, R. P., Pinto, P. A., Rupen, M. P., Smith, 

R. C, Gunn, J. E., & Schneider, D. P. 1991, ApJ, 372, 531 
Li, H. & McCray, R. 1992, ApJ, 387, 309 
Maeda, K., et al. 2008, Science, 319, 1220 

Maeda, K., Nomoto, K., MazzaU, P. A., & Deng, J. 2006, ApJ, 640, 
854 

Maeda, K., ct al. 2007, ApJ, 666, 1069 

Matheson, T., Filippenko, A. V., Ho, L. C, Earth, A. J., & 
Leonard, D. C. 2000, AJ, 120, 1499 



Matheson, T., Filippenko, A. V., Li, W., Leonard, D. C, & Shields, 

J. C. 2001, AJ, 121, 1648 
Matheson, T., ct al. 2008, AJ, 135, 1598 
MazzaU, P. A., ct al. 2007, ApJ, 661, 892 
MazzaU, P. A., ct al. 2005, Science, 308, 1284 
Modjaz, M. 2007, PhD thesis. Harvard University 
Modjaz, M., et al. 2006, ApJ, 645, L21 
Morris, T. & Podsiadlowski, P. 2007, Science, 315, 1103 
Scheck, L., Kifonidis, K., Janka, H.-T., & Miillcr, E. 2006, A&A, 

457, 963 

Schlegel, E. M. & Kirshner, R. P. 1989, AJ, 98, 577 
Schmidt, G. D., Weymann, R. J., & Foltz, C. B. 1989, PASP, 101, 
713 

SoUerman, J., Leibundgut, B., & SpyromiUo, J. 1998, A&A, 337, 
207 

SpyromiUo, J. 1991, MNRAS, 253, 25P 
— . 1994, MNRAS, 266, L61 

SpyromiUo, J., Stathakis, R. A., Cannon, R. D., Waterman, L., 

Couch, W. J., & Dopita, M. A. 1991, MNRAS, 248, 465 
Tominaga, N. et al. 2005, ApJ, 633, L97 
Tuohy, L R. & Dopita, M. A. 1983, ApJ, 268, Lll 
Uomoto, A. & Kirshner, R. P. 1986, ApJ, 308, 685 
Valenti, S., et al. 2008, MNRAS, 383, 1485 
Woosley, S. E. & Bloom, J. S. 2006, ARA&A, 44, 507 



